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Abstract 

We show that the data on 7p and 77 interactions can be derived from the pp 
and pp forward scattering amplitudes using vector meson dominance and the additive 
quark model. The nucleon-nucleon data are parameterized using a model where high 
energy cross sections rise with energy as a consequence of the increasing numbers of 
soft partons populating the colliding particles. We present detailed descriptions of the 
data on the total and elastic cross sections, the ratio of the real to imaginary part of 
the forward scattering amplitude, and on the slope of the differential cross sections for 
pp, pp, 7p, 77, 7p — > jV and 77 — ► V%Vj reactions, where V = p,uj, 4>. We make a wide 
range of predictions for future HERA and LHC experiments and for 77 measurements 
at LEP. 
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1 Introduction 



We show that the data on 7]? and 77 interactions can be derived from the pp and pp forward 
scattering amplitudes using vector meson dominance (VMD) and the additive quark model. 
We first show that the data on the total cross section, the slope parameter B and the ratio of 
the real to imaginary part of the forward scattering amplitude p for pp and pp interactions, 
can be nicely described by a model where high energy cross sections rise as a consequence of 
the increasing numbers of soft partons populating the colliding particles |IJ. The differential 
cross sections for the Tevatron and LHC are predicted. Using this parameterization of the 
hadronic forward amplitudes, we calculate the photoproduction cross sections, slope and p 
value from VMD and the additive quark model. We then obtain 77 cross sections which are 
again parameter- free, demonstrating the approximate validity of the factorization theorem. 

All cross sections will be computed in an eikonal formalism guaranteeing unitarity through- 
out: 



Here, \ is the complex eikonal ( \ = x R +%) j an d b the impact parameter. The even eikonal 
profile function ^ even receives contributions from quark-quark, quark-gluon and gluon-gluon 
interactions, and therefore 



where is the cross sections of the colliding partons, and W(b; fi) is the overlap function 
in impact parameter space, parameterized as the Fourier transform of a dipole form factor. 
This formalism is identical to the one used in "mini-jet" models, as well as in simulation 
programs for minimum-bias hadronic interactions such as PYTHIA and SYBILL. 

In this model hadrons asymptotically evolve into black disks of partons. The rising cross 
section, asymptotically associated with gluon-gluon interactions, is simply parameterized 
by a normalization and energy scale, and two parameters: fi gg which describes the "area" 
occupied by gluons in the colliding hadrons, and J(= 1 + e). Here, J is defined via the 
gluonic structure function of the proton, which is assumed to behave as l/x J for small x. It 
therefore controls the soft gluon content of the proton, and it is meaningful that its value 
(e ~ 0.05) is consistent with the one observed in deep inelastic scattering. The introduction 
of the quark-quark and quark-gluon terms allows us to adequately parameterize the data at 
all energies, since the "size" of quarks and gluons in the proton can be different. We obtain 
fi qq = 0.89 GeV, and fi gg = 0.73 GeV. This indicates that gluons occupy a larger area of the 
proton than do quarks. 

The photoproduction cross sections are then calculated from this parameterization of 
the hadronic forward amplitudes, assuming vector meson dominance and the additive quark 
model. To this end , we introduce -Phad, the probability that the photon interacts as a hadron. 
We will use the value Phad — 1/240 which can be derived from vector meson dominance. Our 
results show that its value is indeed independent of energy. It is, however, uncertain by 30% 
because it depends on whether we relate photoproduction to 7r-nucleon or nucleon-nucleon 
data (in other words, nN and NN elastic cross sections only satisfy the additive quark model 




(1) 



s,b) 



Xqq(s, b) + Xqg(s, b) + Xgg(s, 6) 



i a qq {s)W(b; ji qq ) + a qg (s)W(b; yTW^) + ^ gg (s)W(b; n gg ) 



(2) 
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to this accuracy). Subsequently, following reference 0, we obtain yp cross sections from the 
assumption that, in the spirit of VMD, the photon is a 2 quark state, in contrast with the 
proton which is a 3 quark state. Using the additive quark model and quark counting, the 
7P total cross section is obtained from the even eikonal for pp and pp by the substitutions: 

2 

a ij — * 3 a ij ' 

We will thus produce a parameter-free description of the total photoproduction cross 
section, the phase of the forward scattering amplitude and the forward slope for yp — > Vp, 
where V = p,u,<j). Interestingly, our results on the phase of Vp — > Vp are in complete 
agreement with the values derived from Compton scattering results (7 + p — > 7 + p) using 
dispersion relations. We also calculate the total elastic and differential cross sections for 
■yp — ¥ Vp. This wealth of data is accommodated without discrepancy. 

The 77 cross sections are derived following the same procedure. We now substitute 

4 

&ij * 9 ®ij j 

into the nucleon-nucleon even eikonal, and predict the total cross section and differential 
cross sections for all reactions 77 — > VVj at a variety of energies, where V = p, uj, (f>. 

The high energy 77 total cross section [[J have been measured by two experiments at LEP. 
These measurements yield new information on its high energy behavior at center-of-mass 
energies in excess of yfs = 15 GeV. However, the two data sets unfortunately disagree. We 
here point out that our analysis nicely accommodates the L3 measurements. The analysis is 
sufficiently restrictive to exclude the preliminary OPAL results J|]. It is interesting to note 
that the small eikonal found in our model by fitting n-n data — as shown later — enforces 
naturally the validity of the factorization theorem, 

@pp @yp 

'*yp cy^i^y 

a result independent of the details of our model. We find that the L3 data satisfy the 
factorization theorem whereas the preliminary OPAL data do not. VMD and factorization 
are sufficient to prevent one from adjusting Phad, or any other parameters, to change this 
conclusion. 

An interesting theoretical issue emerges when it is noticed that we applied the additive 
quark model to the full hadronic eikonal, not just to the quark subprocesses in Eq. 2. This 
was not a choice — we found that the data clearly enforced it. For example, if we do not 
apply the quark counting rules to the gluon-gluon subprocess, the model fails to reproduce 
the forward slope of the jp —>■ yV reactions, as well as the ratio of the imaginary to real 
part of the 7p — > jp forward amplitude. It, in fact, fails dramatically. This result may 
suggest either a static structure of the nucleon where the gluons cluster around the original 
valence quarks, as in the valon model 0, or a dynamic picture in which the gluons are 
associated with the interaction of the valence quarks during the hadronic collision. This 
picture is reinforced by correlation measurements between quarks and gluons, derived from 
the observation of multiple parton final state in hadron collisions 0. 
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(3) 



We further emphasize the energy-independence of Phad- in other words, the shapes of 
the total cross sections for ^p and 77 reactions as a function of energy are completely fixed 
by the shape of the nucleon-nucleon data. 



2 High Energy pp and pp Scattering 

In this Section, we will discuss high energy pp and pp scattering. In Section [2.1| , we will discuss 
the theory of our QCD-inspired eikonal and its implementation in fitting the experimental 
data for er tot , P and B to determine the parameters of the model. 



In Section [2.2| , we will compare the experimental data for the elastic scattering cross 
section, elastic as a function of energy with our predictions. 

In Section |2.3| , we will show predictions for differential elastic scattering at yfs = 1800 
GeV, compared to experimental data, and finally, a prediction for the differential elastic 
scattering at the LHC. 

2.1 The Fit to High Energy pp and pp Scattering Data 

We will fit all available high energy forward scattering data above 15 GeV, using bothpp and 
pp, for 

1. (Xtot, the total cross section, 

2. p, the ratio of the real to the imaginary part of the forward scattering amplitude, 

3. B, the logarithmic slope of the differential elastic scattering cross section in the forward 
direction. 

We insist that our QCD-inspired model satisfies: 

1. crossing symmetry, i.e., be either even or odd under the transformation E — > —E, 
where E is the laboratory energy. This allows us to simultaneously describe pp and pp 
scattering. 

2. unitarity. We will use an eikonal formalism to guarantee this. 

3. analyticity. We need this to calculate the phase of the forward scattering amplitude 
and hence, the p value. 

4. the Froissart bound. Asymptotically, we expect that the total cross section will rise as 
(logs) 2 . 

The formalism needed to derive 0tot> P and B from an eikonal are given in Appendix [A|, in 
sections |A.1.1| , |A.1.5| and |A.1.6| , in eqns. (|A8|) , QA14 ) and ( |A20|) , respectively. Details on 



the analyticity are given in Ap pendices p.3.1| and p.4j . The even portion (under crossing) 



of our QCD-inspired eikonal, as noted in eq. (|2]), contains quark-quark, quark-glue and 
glue-glue components, of which the glue-glue portion dominates at high energy. A detailed 
parameterization of this portion is given in Appendix |B.1| , sections |B.1.1| and [B.1.1.2j . The 
even eikonal we finally use is given in eq. (|B23| ). 
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We show in Appendix B.l.l.lj that the Froissart bound is satisfied by our glue-glue 
interaction, and that, asymptotically, the cross section is given by 



e \ . 2 s 



a tot = 2nl — \ log 2 -, (4) 

\figgj S 



as seen in eq. (|B14| ). The parameter e is defined via the gluonic structure of the proton, 



which is assumed to behave as l/x 1+e , for small x. The mass p gg describes the area occupied 
by the gluons in the colliding hadrons. Both are fitted from experiment, with e ~ 0.05 and 
figg m 0.73 GeV. These two parameters, along with the threshold mass m ~ 0.6 GeV and 
the strength of the glue-glue interaction, C gg , are all that is required to specify the glue-glue 
interaction. These 4 parameters dominate the high energy behavior of the nucleon-nucleon 
cross section and are the critical elements of our fit. 

The quark-quark and quark-glue portions are discussed in Appendices |B.2| and [B.3|, and 



are simulated by a constant strength C, a Regge descending trajectory strength CR e g g e even, 
a strength C qg i og for a log term, an energy scale y/s for the log term and a quark size p qq , 
as detailed in eq. ( |B22| ) in Appendix |B.3.1| . We show how to make the even eikonal analytic 



in section |B.3.1| . Details of the odd eikonal are given in Appendix [B.4| , and the odd eikonal, 
which contains two parameters, a strength C dd and a size /i dd, is given in eq. ( |B24j ). 

In all, 11 parameters are used in the theoretical model. The low energy region, for 
-Js <, 25 GeV, where the differences between pp and pp scattering are substantial, largely 
determine the 7 parameters necessary to fit the odd eikonal and the quark-quark and the 
quark-glue portions of the even eikonal. Thus, they largely decouple from the high energy 
behavior, which depends on only 4 of these quantities. Hence, for ^fs > 25 GeV, where there 
is little difference between pp and pp scattering, we really need only 4 parameters for our 
model. 

We use in the fit all of the highest energy cross sections available, (E710§, CDFg and 



the unpublished Tevatron value |10|) which anchor the upper end of our cross section curves. 
The results of the fit are shown in Figs. |l], ^| and || The total cross sections er tot for pp 
(dotted line and circles) and for pp (solid line and squares) are plotted against y/s, the cms 
energy, in Fig. |]. The p values (the ratio of the real to the imaginary portion of the forward 
scattering amplitude) are plotted in Fig. |^ using the same conventions, and the nuclear slope 
B values are similarly plotted in Fig. |3|. 

It can be seen from these figures that we have a quite satisfactory description of all 3 
quantities, for both pp and pp scattering. The \ 2 °f the fit is reasonably good (considering 
the large spreads in the experimental data — in B, in particular, as well as the discrepancies 
in t he highest energy cross sections), giving a \ 2 of 130.3, where 75 was expected. The cross 
section fit of Fig. [l| splits the difference between the values at = 1800 GeV. From Fig. ^ 
we note that the fit to p is anchored at y/s = 550 GeV by the very accurate measurement |p~2f 



of UA4/2 and passes through the E710 point |TT 



The statistical uncertainty of the fitted parameters is such that at 25 GeV, the cross 
section predictions are statistically uncertain to ~ 1.3%, at 500 GeV are uncertain to ~ 1.6% 
and at 2000 GeV are uncertain to ~ 2.5% . 
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Figure 1: The total cross section otot, in m b vs. y/s, in GeV, for pp and pp scattering. The solid line and 
squares are for pp and the dotted line and circles are for pp. 




Figure 2: The ratio of the real to imaginary part of the forward scattering amplitude, p vs. y/s, in GeV, 
for pp and pp scattering. The solid line and squares are for pp and the dotted line and circles are for pp. 
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Figure 3: The nuclear slope parameter B, in (GeV/c) -2 vs. ->/s, in GeV, for elastic pp and pp scattering. 
The solid line and squares are for pp and the dotted line and circles are for pp. 



2.2 Predictions for Elastic Scattering Cross Sections 

We now have all the parameters needed to specify our eikonal. In Fig. f| we have plotted our 
prediction for the elastic cross section <7 e i ast ic, in mb vs. the cms energy \/s, in GeV, along 
with the available data for both pp and pp. The agreement is excellent. 

We note that elastic is rising more sharply with energy than does the total cross section 
a to . Comparing Fig. [I] with Fig. |, we see that the ratio of the elastic to total cross section 
is rising significantly with energy. The ratio is, of course, bounded by the value for the black 
disk[|l4], [0|, i.e., 0.5, as the energy goes to infinity. 

2.3 Predictions for pp Elastic Differential Scattering Cross Sec- 
tions at 1800 GeV 



^From eq. (|A.1.4j) , we can now calculate 4|, the elastic differential cross section as a function 



of for various values of a/s. The calculated differential cross section at the Tevatron 
(y/s = 1800 GeV) is shown in Fig. [5] and compared with the experimental data from E710 [|i~3 |. 
The agreement over 4 decades is striking. 

2.4 Predictions for pp at the LHC 

With the parameters we obtained from our fit, we predict the total cross section at the LHC 
(14 TeV) as er tot = 108.0 ± 3.4 mb, where the error is due to the statistical errors of the 
fitting parameters. 
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Figure 4: Elastic scattering cross sections, Ceiastic, in nib vs. y/s, in GeV, for pp and pp scattering. The 
solid line and squares are for pp and the dotted line and circles are for pp. 




Figure 5: The elastic differential scattering cross section in mb/(GeV/c) 2 vs. \t\, in (GeV/c) 2 , for the 
reaction pp — > pp at yfs — 1800 GeV. The solid curve is the prediction at the Tevatron Collider and the data 
points are from E710. 
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Our prediction for the differential cross section for y/s = 14 TeV, the energy of the LHC, 
is plotted in Fig. ^| It will be a challenge to the LHC to measure this cross section and 




Figure 6: The predicted elastic differential scattering cross section in mb/(GeV/c) 2 vs. \t\, in (GeV/c) 2 , 
for the reaction pp — > pp at y/s = 14 TeV, at the LHC. 

to confirm the predicted structure in \t\. In particular, at small \t\, we predict that the 
curvature parameter C (see Appendix | A . 1 . 7| for details) is negative. For energies much lower 



than 1800 GeV, the observed curvature has been measured as positive. For 1800 GeV, we 
see from Fig. |5| that the curvature parameter C is compatible with being zero. Block and 
CahnJTl, [T5| had pointed out that they expected the curvature to go through zero near the 



Tevatron energy an d that it should become negative thereafter. Basically, the argument is 
that experimentally, the curvature for the then available energies was positive. However, it 
was expected that the scattering asymptotically would approach that of a sharp disk. The 
curvature of a (gray) sharp disk[TJ, [T5|] is always negative, C = — -R 4 /192, where R is the 



radius of the disk. Thus, the curvature had to pass through zero as the energy increased. 
They called 'asymptopia' the energy region (energies much larger than the Tevatron Collider) 
where the scattering approached that of a sharp disk . 



3 7£> Reactions 

In our model, when the photon interacts strongly, it behaves like a two quark system. Taking 
the quark model literally, the eikonal for 7^ scattering is obtained by rewriting the even 
eikonal with the substitutions 

_ 2 
ij 3 ij ' 
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2 1 



(5) 



as 



-CT qq {s)W \ b- J^flqq) + \oqg{s)W lb) J^y/Vqqflgg 



+ l°gg( s ) W ( & ; y ^99 



(6) 



3.1 7J9 Total Cross Section Prediction 

Using vector dominance and the jp eikonal of eq. (H), we can now write, using eq. 



a. 



■yp i 

tot l 



2-Phad / 



-x t "(b,i) cos(x R 7P (&,s)) 



d 2 b, 



(7) 



where Ph a d is the probability that a photon will interact as a hadron. We will use the value 
-Phad = 1/240, which is found by fitting the low energy 7p data. This value is very close to 



that derived from VDM. Using (see Table XXXV, p.393 of ref. |T§) f- = 2.2, f - = 23.6 and 



f 2 

_£. 

4tt 



18.4, we find E 1 



47ra 



1/249, where V = p,u(, 



With all eikonal parameters fixed from our nucleon-nucleon fits and our choice of Phad = 
1/240, we can now calculate cr^(s). Our prediction is given in Fig. [j], where cr^(s) is plotted 




10000 
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Figure 7: The total cross section, a% Q t, in mb vs. y/s, in GeV, for scattering. The solid curve is the 
predicted total cross section. 



against the cms energy. The fit reproduces the rising cross section for jp, using parameters 
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fixed by nucleon-nucleon scattering. We comment here that this prediction only uses the 9 
parameters of the even eikonal, of which but 4 are important in the upper energy region. 
The accuracy of our predictions are pa 1.5%, from the statistical uncertainty in our eikonal 
parameters. 

3.2 'Elastic' Scattering 

We consider as 'elastic' scattering the three vector reactions 

7 + p -> p v i r tual + V -> P + P 

7 + p -> u; virtua i + p -> u + p 

7 + p -> Virtual + P ~> + P, (8) 

where the photon virtually materializes as a vector meson, which then elastically scatters off 
of the proton. The strengths of these reactions is ~ a, the fine-structure constant, times a 
strong interaction cross section. The true elastic cross section is given by Compton scattering 
on the proton, 7 + p — > 7 + p, which we can visualize as 

7 + P -> Pvirtual + P~>p+p->7+P 

7 + p -> Virtual + p^u; + p^7 + p 

7 + P -> ^virtual + P -> + P -> 7 + P, (9) 

is clearly a 2 times a strong interaction cross section, and hence is much smaller than 'elastic' 
scattering of eq. Thus, we justify using eq. (^) to calculate the total cross section which 
we compare with experiment, since only reactions with a photon in the final state (down by 
pa a) are neglected. 



3.2.1 Predictions of p and B for 'Elastic' 7p Reactions 



Using the philosophy of 'elastic' scattering expressed in eq. 
we can immediately write 



Pis) 



Re{i(/1 



d 2 b} 



e iX lp (b,s) 



), and eqns. (|£TJ) and ( P2Cp 



(10) 



and 



! / (l - e ^ v ^) b 2 d 2 b 
~ 2 / (1 _ e *x 7P (M) ^25 



We see from eqns. (|T0|) and (p]) that the predictions for both p and the slope B are free of 
any P^ad factors and hence are independent of normalization — thus being the same for either 
pp, up or 0p final states. 

In Fig. H the value of p from eq. (|I0D is the solid curve, plotted as a function of y/s. 
Damashek and Gilman(l7j calculated the p value for Compton scattering on the proton, using 
dispersion relations, i.e., the true elastic scattering reaction for photon-proton scattering. We 
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Figure 8: The solid curve is the predicted ratio of the real to imaginary part of the forward scattering 
amplitude for the 'elastic' reactions , 7 + p — > Vi + p scattering amplitude, where Vi is p , u>° or tfi° vs. 
y/s, in GeV. The dotted curve is ratio of the real to imaginary part of the forward scattering amplitude for 
Compton scattering |l7|, 7 + p — > 7 + p, found from dispersion relations. It has been slightly displaced from 
the solid curve for clarity in viewing. 

compare this calculation, the dotted line in Fig. [8], with our prediction of p (the solid line). 
The agreement is so close that we had to move the two curves apart so that they may be 
viewed more clearly — this gives us confidence in our approach. 

In Fig. |] the solid curve is our prediction for the slope B against the energy <Js. The 
available experimental data for 'elastic' pp and up final states are also plotted. Again, the 
agreement of theory and experiment is quite satisfactory. 

We note that the predictions of p and B are very critical to our analysis. We have 
assumed that in some manner, the gluons are "attached" to the quarks — when we have a 
two quark system, such as the photon, the factors of 2/3 times a cross section and \J?>/2 
times a p of eq. (§) in the even eikonal of eq. (H) are the same for glue-glue as for quark- 
quark. If we relax this assumption, and only use these factors in the quark, then we get 
sharp disagreement with our predicted p value — being considerably larger than the Compton 
value. This is further exacerbated in the predictions for B, with slopes from 11 (at 5 GeV) to 
16 (GeV/c)~ 2 (at 80 GeV), which are much larger than the experimental values. We stress 
that these conclusions are independent of our choice of Phad- Thus, our model clearly has 
dynamical consequences whi ch will be discussed in detail later. 
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Figure 9: The predicted nuclear slope parameter B, in (GeV/c) -2 vs. y^i, in GeV, for the 'elastic' reaction 
7 + p — > V% + p, where Vi is p , u> or <f> scattering. The solid curve is the prediction. For the reaction 
7 + p — > p° + p, the inverted triangles are the Zeus data, the circles are the HI data and the triangles are 
the low energy data. For the reaction 7 + p — > lj + p, the squares are the Zeus data. 

3.2.2 Predictions of elastic an d % for 'Elastic' 7^ Reactions 

To find the elastic cross sections o"^f stic and differential cross sections da Vp / dt as a function 
of energy, using eq. (|A|3), we write 



^elastic! 5 ) 



pVp 
-'had 



1 -e 



d 2 o, 



(12) 



where P h J^ is the appropriate probability for a photon to turn into V, where V — p, 00 or 
respectively. 

Similarly, the differential scattering cross section is, using eq. (|A13| ), given by 



d<r v * 



s,t) 



3 Vp 
had 



(13) 



(// Jt , J (qb)(l-e^ P ^)d 2 b 

where t = —q 2 . The same factors P ha ^ are used in eq. fll3"D as in eq. (|12[) . 

Since we normalize the experimental data to the elastic cross section found with x lp , and 

1 + f 2 
not to 2(°"ciastic + ^elastic)) we ^ n( ^ that must multiply all -f- by 1.65. Hence, our effective 



coupling s are 
f 2 

^ = 3.6, 

4:71 



f 2 f 2 

JoJeS - 38.9 and = 30.4. 



47T 



47T 



Thus, we define the P ha ^ in eq. (|13f) and eq. fll2|) as 



Ana 



had />2 
J' 



(14) 



(15) 



Veff 
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In Fig. 10, we show our prediction for the 'elastic' reaction 7 + p — > p° +p, where we plot the 
elastic cross section clastic (^) against the cms energy. The solid curve is the predicted cross 
section, the squares are Zeus data, the circles are H 1 data and the inverted triangles the 
low energy data. In Fig. [11], we show our prediction for the 'elastic' reaction 7 + p — > u +p. 




0.00 



10 



100 



1000 



Vs, GeV 



Figure 10: The 'elastic' photoproduction cross section, cr c iastic, in nib vs. y/s, in GeV, for the reaction 
7 + p — » p° +p. The solid curve is the predicted cross section, the squares are Zeus data, the circles are HI 
data and the inverted tri angles the low energy data. 



The solid curve is the predicted cross section, the circles are Zeus data and the squares are 
the low energy data. The agreement of Figs. (|T0| ) and ([0]) is very good, lending further 
support to the model. 

The predicted differential cross sections, dcr/dt, for the 'elastic' reactions 7 +p — > p° + p, 
7 + p — > uj + p and 7 + p — > + p, for diverse energies, are plotted in Figs. [12|, |l| and 
H|, respec tively. The agreement, in absolute normalization and shape, of the predicted 



differential scattering cross sections with the experimental data for all three light mesons for 
all available energies reinforces even more our confidence in our model of jp scattering. 



3.3 How Large is the p, cu and Contribution? 

We sum all of our predictions for the elastic vector interactions for p, to and <p and divide 
this sum by the ratio of cr e i ast ic/o"tot (obtained from x 7P , using eq. and eq. (|A9|) ). We 
call this quantity the total vector meson con tribution. In Fig. |T5|, we compare this to the 
total 7J9 cross section. We find that the fraction of the total cross section for ^p reactions 
that is contributed by the three light vector mesons (p, u and 0) is ~ 0.60. The remaining 
40% could, in the spirit of VMD, be made up of heavier vector meson states, or perhaps 
could be continuum states, or, indeed, a mixture of both. 
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Figure 11: The 'elastic' photoproduction cross section, cr c i a stic, in /ib vs. yfs, in GeV, for the reaction 
7 + p — > uj + p. The solid curve is the predicted cross section, the circles are Zeus data and the squares are 
the low energy data. 




Figure 12: The predicted differential scattering cross section in /xb/(GeV/c) 2 vs. \t\, in (GeV/c) 2 , for 
the 'elastic' reaction 7 + p — * p° + p. The solid curve and the circles (Ballam et al. data) are at y / s= 4.3 
GeV, the dashed curve and triangles (HI data) are at y/s— 55 GeV and the dotted curve and diamonds are 
at y/s= 73 GeV (Zeus data). 
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Figure 13: The predicted differential scattering cross section in /xb/(GeV/c) 2 vs. \t\, in (GeV/c) 2 , for 
the 'elastic' reaction 7 + p — > + p, at ^=80 GeV. The circles are the Zeus data. 
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Figure 14: The predicted differential scattering cross section in /xb/(GeV/c) 2 vs. \t\, in (GeV/c) 2 , for 
the 'elastic' reaction 7 + p — > (f> + p, at -^=70 GeV. The circles are the Zeus data. 
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Figure 15: The dotted curve is the prediction for the total vector meson photoproduction cross section, in 
mb, i.e., the sum of 7+p — > p° +p, 7+p — > u;+p and 7 +p — > (j>+p, divided by the ratio of elastic to total 
cross se ction (see text). For comparison, we also show the solid line and the data which are the predicted 
and experimental total photoproduction cross sections on protons, in mb, vs. yfs in GeV. 



4 77 Interactions 

In this Section, we consider 77 interactions. As we did for 7p interactions, we will take the 



eikonal x 7P ( s ) b) and again multiply every cross section by 2/3 and multiply each fi by y 3/2. 
Thus, we have 



a qq {s)W i^b; -n qq \ + -o qg {s)W ( b; -^H qq H gg 
4 / 3 

+ gVgg(s)W [b] -figg 



(16) 



4.1 77 Total Cross Section Prediction 

Again, using vector dominance and the 77 eikonal of eq. (|T6|), we can now write, using 



eq. 



a 



77 / 

tot l 



2((Phad) 2 / 



-x : 77 (M CO s(x R 7P (M) 



d 2 b, 



(17) 



where, again, Phad — 1/240 is the probability that a photon will interact as a hadron. In 



Fig. [16] we plot our prediction for a^J t (s) as a function of the cms energy and compare it to 
the various sets of experimental data. It is clear that VMD selects the L3 and casts doubt 
on the preliminary OPAL results[H]. 
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Figure 16: The predicted total cross section, er t ot, in nb vs. y/s, in GeV, for 77 scattering. The data 
sources are indicated in the legend. 

4.2 'Elastic' 77 Reactions 

In this Section, we make prediction for 'elastic' 77 reactions, in which both photons turn into 
vector mesons, and then elastically scatter off each other, i.e., 7 + 7 — > Vi + Vj — > V, + Vj, 
where V,, Vj = p, uj, 0. We consider here the 6 reactions 



7 + 7 - 
7 + 7 - 
7 + 7 " 
7 + 7 " 
7 + 7 " 
7 + 7 



P° + P° 
>p° + u 

>p° + 
-> ijj + uj 

-> LU + 

-><!> + <t 



(18) 
(19) 
(20) 
(21) 
(22) 
(23) 



There currently exist no data for such reactions, but hopefully, there will be in the foreseeable 
future — perhaps these predictions will be useful for experimental planning. 



4.2.1 'Elastic' 77 Cross Sections 

To find the total 'elastic' scattering cross sections, we invoke eq. (A9) of Appendix [A. 1.2 
using the eikonal x 77 of eq. fllBD , multiplied by the factors P^ad^} 



?i D hL. *- e v °"el(s) as 



77 ( 
clastic \ 



9 pVi ~p^i 
had had 

pVi p Vj 
had had 



1 -e 



if i ^ j, 
if i = j- 



(24) 
(25) 
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The factor of 2 in eq. (|24|), where there are unlike mesons in the final state, takes into 
account, for example, that either photon in eq. ([OJ) could turn into a p°. In eqns. (pll) and 

(H), the factor P^ d = V- = p, w, 0, where, from eq. O, = 3-6, % = 38.9, 



and 4^ = 30.4. 



We show in Fig. |17] the predicted cross section for reaction fll8|) , in Fig. 18 for reactions 
and ( POD and finally, in Fig. [T9|, for reactions (j2~T|), and ([23]), as a function of ^/s. 
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Figure 17: The predicted cross section in nb vs. y/s, in GeV, for the 'elastic' reaction 7 + 7 — » p° + p°. 



4.2.2 Slope Parameters 5 



Using eq. (|A20|) and the eikonal x 77 , we predict the nuclear slope parameter for the 'elastic' 
reaction of eqns. flT%|)-(|2"3D , as a function of energy. Of course, the slopes are the same for 



all elastic 77 reactions. The predicted slopes B, a s a function of y/s, are shown in Fig. ^0 



4.2.3 Differential Cross Sections 

Using eq. (|A13| ) and x 77 of eq. (|T6|), we can write 



da 



s,t) 



had had 



pV l pVj 
had had 



47T 



Air 



J (qb)(l-e lxll{b ' s) )d 2 b 
J (qb)(l-e^' nib ^)d 2 b 



if i + J, 
if i = j, 



(26) 
(27) 



where t = —q 2 . The factor of 2 in eq. fl26|), where there are unlike mesons in the final state, 
again takes into account, for example, that either photon in eq. ( |P9"D could turn into a p°. 
In eqns. fl2|) and fl27|), the factor P h ^ d = y^ , V- = p, w, 0, where, from eq. flTJ, 



f 2 

4-7T 



9 « /2eff 



38.9, and 



1 4>c« 
4-7T 



30.4. 



In Fig. [H| we show the predicted differential scattering cross section ^ as a function of 
\t\ for the 'elastic' 77 reactions of eqns. ([TJ), (|19|) , and (pTEf), at ^=5 GeV. The solid curve 
is for the reaction 7 + 7 — > p° + p°, the dotted curve for 7 + 7 — > p° + u; and the dashed 
curve for 7 + 7 — > + uj. 



In Fig. p2| we show the predicted differential scattering cross section % for the 'elastic' 



77 reactions of eqns. (|20|) , fl22|) , and (0), at -^5=5 GeV. The solid curve is for the reaction 
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Figure 18: The predicted cross sections in nb vs. s/s, in GeV, for 77 'elastic' reactions. The solid curve is 



for the reaction 7 + 7 



us and the dotted curve for 7 + 7 — > p 
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Figure 19: The predicted cross sections in nb vs. yfs, in GeV, for 77 'elastic' reactions. The solid curve is 
for the reaction 7 + 7 — > oj + uj, the dotted curve for 7 + 7 — > w + and the dashed curve is f or 7 + 7 — > + 0. 
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Figure 20: The predicted nuclear slope parameter B, in (GeV/c) 2 vs. s/s, in GeV, for the 'elastic ' 77 
reactions 7 + 7 — * p° + p°, 7 + 7 — * p° +oj, 7 + 7 — * lu° +uj°, 7 + 7 — » p +</>, 7 + 7 — * + and7+7 — * <fi + cf>. 




Figure 21: The predicted differential scattering cross section in nb/(GeV/c) 2 us. |t|, in (GeV/c) 2 , for 
the 'elastic' 77 reactions, at -^5=5 GeV. The solid curve is for the reaction 7 + 7 — > p° + p°, the d otted 
curve for 7 + 7 — > p + u and the dashed curve for 7 + 7 — > u> + w. 
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7 + 7 — *■ P° + P°i the dotted curve for 7 + 7 — > p° + lu and the dashed curve for 7 + 7 — > uj + cj. 




Figure 22: The predicted differential scattering cross section =|, in nb/(GeV/c) 2 us. |t|, in (GeV/c) 2 , for 
the 'elastic' 77 reactions, at y/s—5 GeV. The solid curve is for the reaction 7 + 7 — > p + tj>, the dot ted curve 
for 7 + 7 — * to + <p and the dashed curve for 7 + 7 — * + <fi. 



In Fig. ^3] we show the predicted differential scattering cross section % for the 'elastic' 
77 reactions of eqns. fll8|) , fll^D, and (|2T|), at ^=20 GeV. The solid curve is for the reaction 
7 + 7 — > p° + p°, the dotted curve for 7 + 7 — > p° + u and the dashed curve for 7 + 7 — > c<j + cj. 



In Fig. we show the predicted differential scattering cross section -J 1 for the 'elastic 



77 reactions of eqns. (p0|) , (|22]) , and (|23|) , at A/i=20 GeV. The solid curve is for the reaction 
7 + 7 — > p° + p°, the dotted curve for 7 + 7 — > p° + to and the dashed curve for 7 + 7 — > uj + cj. 



In Fig. p5| we show the predicted differential scattering cross section ^ for the 'elastic' 



77 reactions of eqns. dig) , (|T9D , and (|2T|) , at ^=70 GeV. The solid curve is for the reaction 
7+7 — > p°+p°, the dotted curve for 7+7 — > p +c<j and the dashed curve for 7+7 — > u; +C(j . 



In Fig. EB] we show the predicted differential scattering cross section ^ for the 'elastic 



77 reactions of eqns. (p0|), (p2]) , and (]23|) , at A/i=20 GeV. The solid curve is for the reaction 
7 + 7 — > p° + p°, the dotted curve for 7 + 7 — > p° + u> and the dashed curve for 7 + 7 — > a; + w. 

It would be most interesting to be able to measure the predicted \t\ structure shown in 
these differential cross section curves. 
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Figure 23: The predicted differential scattering cross section in nb/(GeV/c) 2 vs. \t\, in (GeV/c) 2 , for 
the 'elastic' 77 reactions, at y/s=20 GeV. The solid curve is for the reaction 7 + 7 — > p° + p°, the dotted 
curve for 7 + 7 — > p° + w and the dashed curve for 7 + 7 — > w + u>. 
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Figure 24: The predicted differential scattering cross section %, in nb/(GeV/c) 2 us. |t|, in (GeV/c) 2 , for 
the 'elastic' 77 reactions, at y / s=20 GeV. The solid curve is for the reaction 7 + 7 — > p° + tfi°, the dotted 
curve for 7 + 7 — > lo° + 0° and the dashed curve for 7 + 7 — > 0° + 0°. 
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Figure 25: The predicted differential scattering cross section in nb/(GeV/c) 2 vs. \t\, in (GeV/c) 2 , for 
the 'elastic' 77 reactions, at v / s=70 GeV. The solid curve is for the reaction 7 + 7 — » p° + p°, the dotted 
curve for 7 + 7 — ► p° + lj° and the dashed curve for 7 + 7 — ► a; + cj . 




Figure 26: The predicted differential scattering cross section in nb/(GeV/c) 2 us. |t|, in (GeV/c) 2 , for 
the 'elastic' 77 reactions, at v / i=70 GeV. The solid curve is for the reaction 7 + 7 — > p° + cf>, the dotted 
curve for 7 + 7 — * a; + and the dashed curve for 7 + 7 — ► + (p. 
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5 Summary and Conclusions 



Our conclusions for the total cross section for 7p and 77 are summarized in Fig. |27|. In order 
to scale nucleon-nucleon, ^p and 77 cross sections to a common curve, we have multiplied 
the 7p cross sections by 1/P n ad (=240) and the 77 cross sections by (1/Phad) 2 (=240 2 ). The 
nucleon-nucleon calculation is made using the even eikonal. For clarity, we have not included 
the Opal 77 experimental data. Basically, both the data and our theory approximately satisfy 
factorization, with 

nn— even 7P 

°tot _ "tot /«o\ 

~JfP ~~ -77 ' \ °> 

°tot °tot 

an immediate consequence of the eikonal being small in the energy region considered (up to 
m 2 TeV). The small eikonal we find is consistent with the Tevatron energy not yet being in 
'asymptopia'. 

All data are in agreement with our QCD-inspired eikonal model, which requires that we 
use the even eikonal. When appropriate factors of 2/3 (for quark counting) are introduced, 
and an energy independent factor Phad = 1/240 is introduced, we have a natural explanation 
of 7p interacti ons. Finally, when we again introduce another factor of 2/3, and the factor 
-^had = (24o a ) ' we can ex pl arn the total 77 cross section, agreeing with the L3 data and 
disagreeing the Opal results. We stress that the jp and the 77 experimental data are 
consistent with both Phad being energy independent and the factorization theorem of eq. (]28D . 

We show that VMD, combined with quark counting, fits all available 'elastic' ^p data. 
It even predicts correctly the phase of the forward scattering amplitude for true Compton 
scattering, 7 + p — > 7+p. Our theory allows us to calculate that the three light vector meson 
s, p, uj and 0, account for about 60% of the total 7p cross section. 

Finally, there are also dynamical consequences of our phenomenological model. We see 
that the gluons are carried along with the quarks, since when we link our 2/3 factors exclu- 
sively to the quark composition of our eikonal, we get strong disagreement with the measured 
nuclear slope parameters B for jp 'elastic' interactions, as well as disagreement between the 
p values that we predict for 'elastic' scattering (such as 7 + p — > p° + p) and the dispersion 
relation calculation |T7J for Compton scattering on a proton, 7 + p — > 7 + p. 



- 25 - 




Figure 27: The solid curve is the fitted total cross section for nucleon-nucleon scattering, using the even 
eikonal . The dotted curve is the predicted total cross section for jp scattering multiplied by 1/Phad (=240). 
The circles are the jp data multiplied by 1/Phad- The dashed curve is the predicted total cross section for 
77 scattering multiplied by (1/Phad) 2 (=240*240). The squares are 77 total cross section data multiplied 
by (1/Phad) 2 - 
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A Eikonal Formulation 



In order to ensure unitarity, we utilize an eikonal formalism, evaluating the eikonal in the 
two-dimensional transverse impact parameter space b. We introduce two eikonals, x even {b, s ) 
and x° dd {b,s), even and odd under crossing, respectively (where the proton labora tory 
energy E — > — E), which are both complex and real analytic. In terms of the even and odd 
eikonals, the eikonals we require for pp and pp scattering are given by 

X PP (b,s) = X even (b,s)-x° dd (b,s) 

X PP (b,s) = x ^(b,s) + x odd (b,s). (Al) 

This work largely follows the procedures and conventions used by Block and Cahn [[l4| . 
In terms of the cm. scattering amplitude f c . m ., the cm. differential elastic scattering cross 
section d ^ a (9), the invariant differential elastic scattering distribution % and the total 
cross section cr tot are given by 

>]G (0) = |/ c . m .| 2 , (A2) 



da 



cm. 



dt ~ F l/c - m -' 2 ' (A3) 

a tot = ^Im/ CJn .(0 = O), (A4) 

where k is the cm. system momentum, 9 is the cm. system scattering angle and t = 
—2k 2 (l — cos#) is the invariant four-momentum transfer. Let a(b,s) be the scattering am- 
plitude in impact parameter space. The cm. scattering amplitude ||14|| is given by 



/ c . m .M) = - / e^ b a(b,s)d% (A5) 



7T 



where d 2 b = 2trb db, and q is a two-dimensional vector in impact parameter space b such that 
q 2 = —t. Let the eikonal x(b, s) be complex, such that 

X(b, s) =x R (b,s) + i% l (b,s). ( A6) 

We define our eikonal x{b,s) so that a(b,s), the (complex) scattering amplitude in impact 
parameter space b, is given by 



a(b,s) = 1(1- 



-(l- e -x i (M)+ix R (M) ) . (A7) 



A.l Forward Scattering Parameters and Cross Sections 

We now calculate the forward scattering parameters and various cross sections, using the 
above eikonal formulation in impact parameter space. 
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A. 1.1 Total Cross Section cr tot (s) 

Using the optical theorem, the total cross section cr to t(s) is given by 

, . 47T, 



(It 



k 



■Im/ C . m .(s, 0) 



4 / Im a(b, s) d 2 b 

1 _ e -xj (M) cos ( Xr ( 6j s ) )] ^ ( Ag ) 

where in eq. ( |A8| ) , we used eq. (|A~5|) evaluated in the forward direction (t = 0) and substituted 
eq. (|A7|) to evaluate the total cross section in terms of the eikonal x{b, s). 



A. 1.2 Elastic Scattering Cross Section cr e i(s) 



^From eq. (|A3|) , we can evaluate the elastic scattering cross section a e \(s) as 

0"elastic(s) = ^ J I /cm. (s,t) | 2 dt 



P / |/am.M)W 

^- 6 \(6, s)a*(6', s) d 2 qd 2 bd 2 V 



(27T) 



7T 



a(b,s)a*(b',s)5 2 (b-b')d 2 qd 2 b 



4 / |a(M)r<™ 



1 _ e -Xi(M)+*X R (M) 



A. 1.3 Inelastic Cross Section <7i n (s) 



Thus, using eq. ( |A~8|) and eq. (|A9|) , the inelastic cross section, defined as <7 tot (s) 
is given by 



^inelastic ) 



,-2Xj(M 



} rf 2 6. 



(A9) 

^elastic (^0 j 

(A10) 



A. 1.4 Differential Elastic Scattering Cross Section %?(s,t) 



A convenient way of calculating the differential elastic scattering cross section of eq. (|A.1.4| ) 
is to note that an alternative way to write eq. ( |A5| ) is to introduce an integral representation 



of J (see eq. (9.1.18), ref. g9 



1 r 2n 

Uz) = — e i2COS ^ 
Ztt Jo 



We can then rewrite eq. (|A5|) as 

poo U 

fc m ( s i t) — 2k / J (qb) a(b, s)b db = — I J (qb) a(b, s) d 2 b. 
Jo IX 



(All) 



(A12) 
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Finally, using eq. ( A12 ) and eq. 

as 



we now write the differential scattering cross section 



da 



J (qb)(l-e- x ^ b ' s)+ix ^ b ' s) )d 2 b 



a more convenient computational form. 



(A13) 



A. 1.5 p(s) 

To calculate p, the ratio of the real to the imaginary portion of the forward nuclear scattering 
amplitude, we write 

Ref c . m .(s,0) 



P(s) 



Im/c. m .O,0) 

Re {i(J 1 - e-*i( 6 > s ) +i *R^) d 2 b} 



Im{z(f(l 



d?b) 



(A14) 



A. 1.6 Nuclear Slope Parameter B 

The nuclear slope parameter B is defined as 

Bis) 



d 
dt 



, da , . 
In-(M) 



(A15) 



t=0 



(A16) 
(A17) 

, we can eventually write the general 

(A18) 



Beginning with eq. (|A5| ) , 

/c. m .M)oc J e^a(b,s) d% 

we expand the exponential about q = to get 

k(s, t)<x J[l + iq-b- ~{q- b) 2 + -- -}a(b, s) d 2 b 

With this expansion and the definition of B in eq. 
expression for B as 

_ Re{/ °° dbba{b,s)f™ dbb 3 a*{b,s)} 

2|/ °° dbba{b,s)\ 2 

If the phase of a (b, s) is independent of b (this is the case when we either have a factoriz- 
able eikonal or an eikonal with a constant phase), eq. ( |A18|) reduces to the more tractable 
expression 

B = SS°db^a(b,s) = J b 2 a(b, s) d 2 b 
2J™dbba(b,s) 2ja{b,s)d 2 b 

We note from eq. (|A19|) that B measures the size of the proton, i.e., B is one-half the average 
value of the square of the impact parameter b, weighted by a(b,s). Again, introducing the 
eikonal into eq. ( |A19|) , we find 



B 



which we use to compute B. 



(A20) 
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A. 1.7 The Curvature Parameter C 

The t dependence of the elastic differential cross section is described at small \t\ as 



da 
~dJ [ 



s,t) 



da 

~dl [ 



s,t) 



da , 



t=0 



t=0 



Bt+Ct 2 



1 + Bt + ( T + c ) f2 + ' 



(A21) 



We state, without proof, that the curvature parameter C is given |]m|| by 



C(s) 



32 



J ^1 _ g-XjCM+iXaCM) d 2 6 



1 

16 



/ (l - e" x i (f, ' s)+ix H (6 ' s) ) 6 2 d 2 b 



/(1 



_ e -X I (fc I s)+«X R (6,s) 



) d 2 b 



(A22) 



where it was again assumed that the phase of the eikonal is independent of b (see eq. (4.45) 
of ref. [14]]). From eq. ( A22|) , we see that the curvature C can be positive, negative or zero, 
whereas the nuclear slope parameter, from eq. ( |A20| ), must be positive. 



B QCD-Inspired Eikonal 

B.l Even Eikonal 

The even QCD-Inspired eikonal Xeven is given by the sum of three contributions, glue-glue, 
quark-glue and quark-quark, which are individually factorizable into a product of a cross 
section a(s) times an impact parameter space distribution function W{b i.e.,: 

x even {s,b) = XggM) + X qg {s,b) + Xqq(s,&) 

= i [a gg (s)W(b ; /i gg ) + a qg {s)W{b ; ^n m n gg ) + a m (s)W(b ; // qq )] , (Bl) 

where the impact parameter space distribution function 

W(b^) = -^-^bfK^b) (B2) 

is normalized so that 

W(b;n)d 2 b= 1. (B3) 



Hence, the ex's in eq. ( pl|) have the dimensions of a cross section. 

The factor i is inserted in eq. ( |B~1| ) since the high energy eikonal is largely imaginary (the 
p value for nucleon-nucleon scattering is rather small). 

As a consequence of both factorization and the normalization chosen for the W(b; n), it 
should be noted that 

/ x ^( s ,b)d 2 b = t[a gg (s) + a qg (s) + a m (s)], (B4) 
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so that, using eq. ( |X8| ) for small x, 

atr(s) =2lm{t [a gg (s) + a qg (s) + a qq (s)]} . (B5) 

In eq. ([BID, the inverse sizes (in impact parameter space) /x gg and /i gg are to be fit by 
experiment, whereas the quark-gluon inverse size is taken as y/fi qq fi gg . 

B.l.l The (T gg Contribution 

Modeling the glue-glue interaction after QCD, we write the cross section cx gg (s) in eq. (]Bl| ) 

as 

v gg (s) = C gg N 2 g J Z gg 6(s - ml)F gg Lx 2 = ^ d (-) , (B6) 



where 



o 



The normalization constant C gg and the threshold m are to be fitted by experiment (in 
practice, the threshold is taken as m = 0.6 GeV and the strong coupling constant a s is 



fixed at 0.5). The constant N g in eq. is given by N g = | (5-e)(4- e )(3- e )(2-e)(i-e) j Using 
the gluon structure function 



we can now write the function F gg in eq. (|B6| ) as 



Fgg = J J fs( x i)fs( x ^( x i x 2 = T ) dxi dx 2 - (B9) 

After carrying out the integrations, we can now explicitly express c gg (s) as a function of s. 
The parameter e in eq. (|B8[) is to be fitted by experiment (in practice, we fix it at 0.05). 



B. 1.1.1 High Energy Behavior of <r gg (s) — the Froissart Bound 

We note that the high energy behavior of c gg (s) is controlled by 
lim / drFaJr) ~ / dr- 



- : - lml/s J mils ' T 1+e 



~\M^u)' (B10) 

where e > 0. The cut-off impact parameter b c is given by 

cW gg {b c]f i gg )s e log{s) ~ 1, (Bll) 
where c is a constant. For large values of fib, we can now write eq. (|B 1 1|) as 

c'G^A) 372 ^ 99 V log(s) ~ 1 (B12) 
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with c' another constant, and therefore, 

b c = — log - + O flog log 4) • (B13) 
Vgg so V s oJ 
We reproduce the Froissart bound from QCD arguments, 

atot = 2tt ( — ) log 2 -, (B14) 

as we go to very high energies, as long as e > 0. The usual Froissart bound coefficient of the 
log 2 j- term, 1/m 2 = 20 mb, is now replaced by (e//i 9S ) 2 ~ 0.002 mb. Note that fi gg controls 
the size of the area occupied by the gluons inside the nucleon. 



B.l.1.2 Evaluation of the a gg Contribution 



In the following, we set the matrices a(0) = — a(5) = —411/10, a(l) = — a(4) = —975/2, 
a(2) = -a(3) = -600 and 6(0) = 6(5) = -9, 6(1) = 6(4) = -225, 6(2) = 6(3) = -900. The 
result is 



a gg (s) = C gg Z gg N 2 g f F gg dr 

J TO 



CggSgg./Vg x 

5 r a u\ _ m . (a{i) - b M bit) 



e e 



-975 1 225 / -975 1 225 00 r- , / \- 

-2" + ~e 1-e — + ~e 225 lo g( r o) 



_|_ £ _ T , _ _ , 

1-e V !- e !- e >/ 

-600 + fg 2-«/ r -gM+jg 9001og(r ) 
+ 2-e T ° v 2-e 2-6 

600 + E _ »-« /600+jg _ 9001og(r ) \ 
+ 3-e V 3 - e 3 " e J 

975 , 225 / 975 , 225 00 r- i / \\ 

T + ~e 4-e — + —e 22 5 log(r ) ^ 

+ — ; — ; r o 



4-e 4-e / 

411 , 9 / 411 , 9 „, / \\l 2 

+ ^±s_^(^±s_£^0)}, whereT0 ^. (B15) 

We note that we must fit the following 3 constants in order to specify a gg . 

1. the normalization constant C gg . 

2. the threshold mass m . 

3. e, the parameter in the gluon structure function which determines the behavior at low 
x (oc l/x 1+e ). 
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B.2 The cr qq Contribution 

If we use the toy structure function 
(1 — x) 3 

AO*) = iL #- > (B16) 



we can write 

f \ m ° 1 S , T> ( m °' 

a m {s) oc —= log — + \—= 



~ constant + (B17) 

where P is a polynomial in mo/ 

Thus, we approximate the quark-quark term by 

Cqq( s ) = S g g + Cpj,egg e even" J=^J j (B18) 

where C and CRegge even are constants. Thus, cr qq (s) simulates quark-quark interactions with 
a constant cross section plus a Regge-even falling cross section. 
We must fit the following 2 constants in order to specify cr qq : 

1. the normalization constant C. 

2. the normalization constant CR eg ge even- 

B.3 The dqg Contribution 

If we use the toy structure function 

f g (x) = 2^2, (B19) 
and the toy structure function f q (x) of eq. ( P16|) we can write 

a qg (s) oc C"log- + C'V'(^) 

« C'log — + C", (B20) 
so 

where C" and C" are constants and T 3 ' is a polynomial in itlq/^/s. 

Thus, if we absorb the constant piece C into the quark-quark term, we can approximate 
the quark-gluon term by 

g 

Oqg (S) = SggC qg log log — , (B2 1 ) 

so 

where C qg \ og is a constant. Hence, we attempt to simulate diffraction with the logarithmic 
term cr qg (s). 



We must fit the following 2 constants in order to specify <7 qg : 

1. the normalization constant C qg i og . 

2. s , the square of the energy scale in the log term of eq. ( B21|) . 
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B.3.1 Making the Even Contribution Analytic 

The total even contribution, which is not yet analytic, can be written as the sum of equations 



BT5), piq and P21j i.e., 



Xeven = i { (T gg (s)W (b ] fi gg ) 



+ S g g ( C + Cftegge even — ^= J W(b ; fiqq) 



+ Sg g C qglog log —W{b; v 7 ^^))- (B22) 



For large s, the even amplitude in eq. (fB22|) can be made analytic by the substitution (see 
the table on p. 580 of reference []I4 |, along with reference |0 



Thus, we finally rewrite the even contribution of eq. ( |B22| ), which is now analytic, as 

Xeven = % j ^gg (se~^ /2 ) W (6 J /i gg ) 

+ S gg + CR e gg e even y=e i7r//4 ^ W(b J // qq ) 

+ SggC qglog (log ^ - W(6 ; VW^)}. (B23) 

To determine the impact parameter profiles in b space, we also must fit the mass param- 
eters /igg and /x qq to the data. We find masses fi gg ~ 0.73 GeV and // qq ~ 0.89 GeV. 

B.4 The Odd Eikonal 

It can be shown that a high energy analytic odd amplitude (for its structure in s, see eq. 



(5.5b) of reference Ulif , with a = 0.5) that fits the data is given by 
X f d (b,s) = ~a odd W(b ]f i odd ) 



-C odd £ g g^e i7r / 4 W(b; /i odd ), (B24) 
v s 



with 



W{b, /i odd ) = ^(/i odd 6) 3 K 3 (fi odd b), (B25) 
normalized so that 

W(b;{i)d 2 b= 1. (B26) 



Hence, the cx odd in eq. ( |B24j ) has the dimensions of a cross section. 

In order for C odd to be positive, a minus sign has been inserted in eq. ( B24|) . 
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With the normalization (eq. ( |B25|) and eq. (|B26|) ) chosen for W(b, /Vid)> we see that 
X odd (s,b)d 2 b = a odd (s), (B27) 



so that, using eq. 



for small x, 



odd/ 



0"t"oTl s J = 2Ima odd (s). 



(B28) 



In order to determine the cross section cx dd, we must fit the normalization constant C Q dd- 
To determine the impact parameter profile in b space, we also must fit the mass parameter 
//odd to the data. We find a mass /x Q dd ~ 0.53 GeV. 



We again reiterate that the odd eikonal, which we see (from eq. (|B24j )) vanishes like 
accounts for the difference between pp and pp. Thus, at high energies, the odd term vanishes, 
and we can neglect the difference between pp and pp interactions . 
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